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Abstract. We study the Coulomb-to-dipole transition which occurs when the 

separation d of an electron-hole bilayer system is varied with respect to the 
characteristic in-layer distances. An analysis of the classical ground state 
configurations for harmonically confined clusters with N < 30 reveals that the 
energetically most favorable state can difi^er from that of two-dimensional pure dipole 
or Coulomb systems. Performing a normal mode analysis for the = 19 cluster 
it is found that the lowest mode frequencies exhibit drastic changes when d is varied. 
Furthermore, we present quantum-mechanical ground states for A'' = 6, 10 and 12 spin- 
polarized electrons and holes. We compute the single-particle energies and orbitals in 
self-consistent Hartree-Fock approximation over a broad range of layer separations and 
coupHng strengths between the limits of the ideal Fermi gas and the Wigner crystal. 
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1. Introduction 

Self-organized structure formation, in particular Coulomb crystallization [1], is among 
the most exciting cooperative phenomena in the field of charged many-particle systems. 
In the case of finite, parabolically confined systems extensive experimental and 
theoretical work on various types of two- and three-dimensional systems has revealed 
that in the strong coupling limit charged particles can arrange themselves in a highly 
ordered crystalline state with a nested shell structure. Examples are ions in Paul 
and Penning traps [2, 3], dusty plasmas [4]-[ll] and electrons in quantum dots and 
wells [12]-[18]. For these also called artificial atoms Mendeleev-type periodic tables 
were found including characteristic occupation numbers, shell closures and unusually 
stable magic configurations. For a recent overview see [1]. Recently there is growing 
interest in two-dimensional (2D) dipolar macroscopic systems [19]-[24] as well as finite 
size dipolar (quantum) clusters in small-scale confinement potentials [25]-[32]. While in 
particular the ground state and dynamical properties of 2D mesoscopic pure Coulomb 
and pure dipole interacting particle ensembles in parabolic confinement potentials are 
well understood, the behaviour of real three-dimensional electron-hole double layer 
systems, where the dipole approximation is not valid, is still poorly investigated. This 
despite the fact that the additional degree of freedom, i.e. the layer separation d, is 
expected to allow for a variety of interesting new efi'ects which are due to the possibility 
of tuning the effective in-layer interaction potential. 

The results presented in this paper arc applicable to semiconductor heterostructures 
and coupled quantum dots as well as to molecular systems, where the dipole moment 
of the charge carriers and thus the interaction strength is tunable, e.g. [32, 33] .J 
For a consistent formulation we concentrate on the problem of two vertically coupled 
symmetric layers containing parabolically confined, spin-polarized electrons and holes 
of identical particle number A^e = Nh = N and effective masses m* — rn^ — m*, 
respectively. The underlying Hamiltonian is 

H ^ He + Hh - He-h , (1) 
with the intra- and interlayer contributions 



^ ^ 47r£y(r, - r,-)2 + d^ 

where the electrons (e) and holes (h) are confined to planes of zero thickness which 
are a distance d apart. The 2D vectors ri(j) are the in-plane projections of the particle 

I Another natural source of confinement arises in low-dimensional semiconductor structures from 
defects and well width fluctuations. This leads to local potential minima for the charge carriers causing 
localization of free and bound charges (excitons, biexcitons, trions), e.g. [34, 35, 36]. 



Coulomb-dipole transition in mesoscopic electron-hole bilayers 



3 



coordinates, e the elementary charge and e the static permittivity. The strength of the 
confinement is controllable by the trap frequency Uq. 

The most fascinating property of this system is that the effective in-layer particle 
interaction changes with the interlayer separation d : from Coulomb interaction at 
large d, where both layers are decoupled, to dipole interaction at small o? — > 0, where 
the attractive interlayer interaction leads at low temperature to vertical electron-hole 
coupling and formation of vertically aligned dipoles — excitons. On the other hand, at 
intermediate values of d, when the repulsive intra- and attractive interlayer interaction 
energies according to equations (2) and (3) are comparable, the system shows a real 
three-dimensional behaviour. In reference [28] it was reported that, as a consequence of 
the Coulomb-dipole transition, the considered system can exhibit structural changes of 
its ground state shell configuration when d is varied. 

In section 2 we extend these results and present a systematic study of the classical 
ground states, varying d for mesoscopic clusters with iV < 30 particles in each layer. 
Further, we extract the fundamental dynamical features in the case of weak excitation 
by solving the dynamical (Hessian) matrix for the ground state configurations found 
in section 2. Doing this, in section 3 we discuss the d-dependence of the collective N- 
particle modes for the = 19 cluster. Here we highlight the close relationship between 
structural and collective dynamical cluster properties as rotation of shells and vortices. 
In section 4 we extend the analysis to fcrmionic e-h quantum bilayers utilizing a self- 
consistent Hartree-Fock ansatz. In particular, Coulomb-to-dipole transition induced 
(critical) quantum phenomena arc presented for the clusters with = 6, 10 and 12 
electrons and holes. The results include the A^-particle densities and the single-particle 
spectrum and orbitals as function of coupling strength A and layer separation d. 



2. Classical ground state transitions 

The classical ground state corresponding to the equations (1) to (3) is described by the 
Hamiltonian H = He + Hh — H^-h without the kinetic energy, i.e. 

This dimensionless form is obtained applying the transformation rules {r — > r/ro, 
> E/Eq, d— > d/ro} with the characteristic length tq = (e^/27r£ma;Q)^/^ and energy 
Eq = (ma;oe^/327r^£^)^/^. Note that model (4) contains no exphcit dependence on 
the trap frequency uoq. The considered classical model system in the ground state 
is completely defined by only two parameters: the particle number N and the layer 
separation d^ which also infiucnces the in-layer density. 

The ground state configuration is the energetically lowest of all possible stable 
states, whose number rapidly increases with N , and all these have to be found and 
checked. This task is complicated, since many of the different stable states are 
energetically close, requiring high-accuracy computations. A systematic search for 
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the global minimum-energy structure in the 4A'^-dimensional configuration space was 
performed by means of an optimized molecular dynamics annealing technique utilizing 
an adaptive step size control [8, 28]. For each value of and d the annealing process 
was repeated for a large (A^- and rf-dependent) number of times. At this a slow (long) 
annealing process ensures to find the lowest-energy state with high probability. The 
critical points of structural transitions dcr were identified as crossing points of the 
energies of the lowest-energy states as function of layer separation d. 

Extending the analysis of [28] we obtained a periodic table for the particle numbers 
A^ < 30 including all structural transitions occurring when d is changed, see table 1. 
In the limits of pure dipolc and Coulomb interaction our results are in full agreement 
with those of reference [25] and references [12, 16], respectively. § Analyzing the clusters 
A^ < 18 only transitions for A^ = 10 and A^ = 12, reported in reference [28], are found. 
Due to the much larger configurational space, and thus accordingly higher number of 
low-energy metastable states, for the clusters A^ = 19 ... 30 in total 6 particle numbers 
reveal Coulomb-dipole transitions: N — 19, 21, 23, 26, 29 and 30. In particular, two 
transition types are identified: 

(A) While for the majority of the investigated clusters the ground state shell 
configuration of the single layer Coulomb and dipole case are identical, for A^" = 10, 
21, 23, 26 and 29 this is not the case. When changing from a long-range Coulomb 
to a short-range dipole interaction a higher particle number on the inner shell 
becomes favourable. A similar trend is also known from 2D [5, 6] and 3D [9, 10] 
Yukawa-clusters when the screening strength is increased. || 

(B) A second type of transition is found for A^ = 12, 19 and 30 that cannot be concluded 
from diff'erent shell occupations in both limits of ci: At large values of d again a 
transition of type (A) takes place, which increases (decreases) the particle number 
on the inner (outer) shell when d is reduced. But interestingly, at small values 
of d a second kind of transition to a sixfold-coordinated, commensurate particle 
configuration is found allowing for an energetically more favourable closed packing 
of the composite dipoles. Such symmetry-induced re-entrant configuration changes 
are only observed in cases where highly symmetric, "magic" configurations with a 
bulk-like triangular structure are involved. 

These findings coincide with those for single layer statically screened Coulomb systems. 
Here a change from the long-range Coulomb towards a short-range Yukawa potential 
by variation of the screening length leads to analogue ground state transitions for the 
particle numbers A^ = 10, 12, 19 and A^ = 21, 23, 26, 29 as reported in reference [6] and 
[5] , respectively. Further, a comparison of the ground and metastable states of the single 
layer Coulomb system (cf. table 1 in reference [16] for A^ < 30) shows that if and only if 

§ In reference [12] the ground state for = 29 was erroneously given as (5,10,14). This was corrected 
in reference [16]. 

II The effect is due to the radial balance of total internal and external F"^^* forces on each particle. 
In contrast to Coulomb, short-range (dipole or Yukawa) forces do contribute to F"^^* which requires a 
higher density towards the center to stabilize the cluster matching F™* = F^^*. For details see [10]. 
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Table 1. Ground state shell structures for 2D Coulomb, bilayer and dipole clusters 
of N particles in a parabolic confinement. The arrows indicate the direction of the 
ground state transition from large to small values of d. Magic (commensurate) shell 
configurations are underlined. For N < 5 only a single shell is populated for all values 
of d. For all configurational transitions the critical layer separation dcr as well as 
the corresponding total energy per composite dipole E^r/N is given. Note that the 
binding energy 1 /d which ensures the exact vertical alignment of the electron-hole pairs 
is excluded from the energy values as it is independent of the cluster configuration. 



an energetically close metastable configuration with higher center particle number than 
the ground exists, in fact, a transition of type (A) in the corresponding bilayer system 
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is found. This underlines the Coulomb-to-dipole transition induced density change 
effecting configurational transitions of type (A). In contrast, transitions of type (B) are 
geometry induced supporting an equally distant, closed packed particle arrangement. 

Among all transitions, the most interesting arc those of type (B). As an example, 
we study the N = 19 chister. Here between d = 0.417 and d = 2.182, the "magic" 
configuration (1,6,12) is replaced by the configuration (1,7,11) which possesses a much 
lower orientational order [37]. Therefore, it is interesting to analyze the normal modes 
of this cluster and their dependence on d. 

3. Collective A^^-peirticle modes 

Starting from the ground state configurations given in section 2, we are interested in the 
collective excitation behaviour in dependence on d. Here we will focus on the cluster 
with = 19 where, upon changing d, finite size effects are expected to play a key 
role as the ground state structure changes between the hexagonally ordered (1,6,12) 
configuration and the (1,7,11) circular ring structure as discussed in section 2. 

To derive the dynamical properties in the limit of weak excitations we perform a 
normal mode analysis [27, 37, 38, 39]. For small particle displacements u(t) = r{t) — R 
around their ground state position R, expansion of the potential energy U, equation 
(1), around R leads to 



where Uq is the minimum potential energy. In the stationary states the linear (force) 
term vanishes and the second-order partial derivatives provide the elements Hij of the 
2 X 2A^ Hessian matrix. In the frame of the harmonic approximation the resulting cluster 
dynamics is given as a superposition of these collective (normal) modes statistically 
weighted according to the eigenvalues of Ti. which are proportional to the squared mode 
oscillation frequencies uf. In the following these eigenfrequcncies will be given in units 



3.1. Classification of normal modes 

As a result of the eigenmode computation we obtain for each stable configuration of 
the A" = 19 cluster a complete set of 76 eigenvalues and eigenvectors. A selection 
of characteristic and energetically low-lying eigenvectors for d — dcr — 2.182, i.e. 
intermediate between Coulomb and dipole regime, is given in figure 1. As shown in 
reference [27], in dipolar bilayer systems the total number of modes can be divided in 
two types which will be distinguished by the following nomenclature: 

("•") labels modes with in-phase collective particle motion in both layers, see figure 1 a), 
d), e), h) to o), and 






Coulomb- dipole transition in mesoscopic electron-hole bilayers 
Selected normal modes of the (1,7,11) configuration 



(a) < 5 • 10- 

(f) = 0.4548 



(b) ^_ = 2.2 • 10-* (c) = 0.0048 (d) uj' 





(g) Lol^_ = 0.4549 (h) c^^^^ 





r 




4^ ^4—*' V*. 



Selected normal modes of the (1,6,12) configuration 



(k) tj^2+ = 0.5866 

t ^ * M 



(1) = 0.6507 

• T * 



V2+ 



: 0.4519 



0.9067 (i) a;|+ = 2.000 



0.4520 



(j) = 6.8703 



1.3819 (n) ^2^+ 



6.8704 



t ^ • i . ^ 



1.9842 (o) 



Figure 1. Top view of tlie eigenvectors of selected characteristic and low-energetic 
normal modes for the = 19 cluster at d = dcr — 2.182 (ordered by frequency, cf. 
numbers above the figures). The points mark the particle positions. The different 
shaped (and colored) arrow heads are assigned to the normal mode eigenvectors in the 
two different layers and indicate direction and amplitude of particle motion. Modes 
with in-/ anti-phase motion of 2 layers are labeled with a +/~ sign, respectively . 
Top rows: Eigenvectors of the (1,7,11) configuration: a) inter-shell rotation (SR"*"), b) 
anti-phase layer rotation (LR~), c) anti-phase inter-shell rotation (SR~), d) and e) 
in-phase vortex pairs (V2+), f) and g) anti- phase vortex pairs (V2^), h) asymmetric 
in-phase 4- vortex mode (^4+), i) sloshing mode {S^), j) breathing mode (S^). 
Bottom row: Eigenvectors of the (1,6,12) configuration: k) in-phase vortex pair (V2+), 
1) in-phase inter-shell rotation (SR+), m) in-phase 4- vortex mode (F4+), n) in-phase 
transverse surface wave (SW^), o) breathing mode (B^). 



( ) labels modes with anti-phase motion of both layers, see figure 1 b), c), f) and g). 

Consider first the top rows of figure 1 which show the eigenvectors of the normal 
modes of the (1,7,11) configuration. The energetically lowest collective particle motion 
is in all cases the center of mass cluster rotation mode — the in-phase layer rotation 
LR"*". The eigenfrequency of this directed rotation is u; = as for this motion there is no 
restoring force. Beside this (trivial) mode there are three additional rotational modes: 
a) inner versus outer inter-shell rotation SR^, b) the anti-phase rotation of both layers 
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LR~ and c) anti-phase inter-shell rotation SR". 

Another set of low frequency modes are four vortex pair modes: d) in-phase vortex 
pair V2"'" and e) perpendicular oriented vortex pair V2^, f ) and g) two anti-phase vortex 
pairs. In the present, isotropically confined 2D system rotationally asymmetric modes 
are typically two-fold degenerate with respect to the spatial alignment of the vectors, cf. 
d), e) and f), g), respectively. This leads to the fact that, taking into account the two 
possible phasings of relative particle motion in both layers, a majority of mode types 
occures as a set of four. Considering this, in the following only one mode per set of four 
is shown as for the rotational asymmetric, low-energy mode h) which has the interesting 
feature that it supports a single particle exchange between the inner and outer shell, i.e. 
a transition from the (1,7,11) to the (1,6,12) configuration. 

In the case of pure radial eigenvectors, such as the (in-phase) breathing mode j) as 
coherent radial motion (compression/expansion) of all particles, there exists one pair of 
modes only. In addition to j) there is an anti-phase breathing mode B~ with frequency 
a;|_ = 7.9522. Further "universal modes" that are independent of particle number and 
configuration is the center of mass sloshing mode i) with trap frequency ujo which has a 
corresponding anti-phase shear or dipole oscillation mode S~ (each two-fold degenerate). 

For all these modes a corresponding mode of the (1,6,12) configuration is found. 
In particular: k) the V2"'"-mode, 1) the mode of inter-shcll rotation SR"^, m) an 
energetically low V4+-mode, here supporting a center directed transition of a particle 
on the outer shell, and further two examples of radial modes, n) a transverse surface 
wave and o) the breathing mode. 

3.2. Change of normal mode spectrum with layer separation 

After the classification of the collective modes we now consider the oscillation frequency 
dependence on the layer separation d of the = 19 cluster, sec figure 2. Of special 
interest are thereby the two configuration changes of the ground state and their effect 
on the collective dynamical cluster properties. 

Starting at small values of d, an increase of the e-h separation leads to a growing 
cluster size due to a stronger in-layer particle repulsion resulting from a change of the 
effective interaction from dipole to Coulomb. This implicates a gradual decrease of the 
mode eigenfrequencies with d since the coupling of all 2A^ particles becomes less rigid 
and the restoring forces weaken. Only the two-fold degenerate center of mass oscillations 
are found to be constant at c<j|+ = 2, independent of the interlayer coupling strength or 
even configuration changes. Confirming reference [27], the breathing frequency gradually 
proceeds from a;^.,. = 10 in the limit of dipoles (rf ^ 0), to a value of 0;^+ = 6 in the 
limit of decoupled layers {d — > 00). Moreover, modes supporting a transition from the 
(1,6,12) to the (1,7,11) state and vice versa, i.e. the eigenmodes h) and m) in figure 1, 
are found at low frequencies, i.e., at low excitation energies. 

As discussed in section 2, the ground state transitions for the A^ = 19 cluster 
occur at the critical values of d^-^ = 0.417 and dir'* = 2.182 and are accompanied 




Figure 2. Complete normal mode spectrum for = 19 as function of layer separation 
d. At dir^ — 0.417 the ground state configuration changes from (1,6, 12) to (1,7, 11) 

(2) 

and at diV = 2.182 from (1,7, 11) to (1,6, 12) resulting in a qualitative change of the 
mode frequencies. The eigenvectors of the selected modes are visualized in figure 1. 
Modes with in-phase (anti-phase) oscillation of both layers are plotted with dashed 
(solid) lines. Note that the SR+ mode continues in the range 1 < d < 2.182 with a 
value smaller than 5 • 10^^. For notation of modes see figure 1. 
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by abrupt spectrum transformations. The strongest effect is observed for the in- 
phase inter-shell rotation SR+ with a remarkable jump of the mode frequency ^|^+ 
by more than four orders of magnitude. This decrease can be explained by comparing 
the SR'^ mode eigenvectors of the (1,7,11) and (1,6,12) configurations, see figure 1 a) 
and 1). In the latter case, the oscillation vectors of all particles on the inner shell 
are directed towards particle positions on the outer shell which strongly increases the 
restoring forces in the case of the (1,6,12) configuration resulting in a much higher 
frequency u}^^^+^^ than ijU^g^+^\ The exceptional low frequency 00^^^+^^ agrees with 
results for single layer Coulomb crystals. In reference [37] the minimal (non-zero) 
excitation frequency for (1,7,11) and the comparable, non-magic (1,7,12) configuration 
was reported to be that of the inter-shell rotation with Ugj^ ~ 10~^. Confirming this, 
quantum Monte Carlo simulations [14] revealed that the orientational inter-shell melting 
temperature of the incommensurate (1,7,12) configuration is much lower than for the 
high symmetric (1,6,12) structure. In particular, a 9 (!) orders of magnitude difference 
of the orientational melting temperatures and critical densities of both configurations 
was found. This shows that the given classical results are of practical relevance also for 
quantum systems at moderate densities. 

Moreover, with respect to the dipole-to-Coulomb transition we found that in 
the dipole regime at small d the corresponding modes with in-phase and anti-phase 
oscillation of both layers are energetically clearly separated, cf. SR^ and V2^ in figure 
2. Energetically lowest arc the two (degenerate) in-phase vortex pair oscillations V2+. 
With a gradual transition to the limit of uncoupled layers, the c-h attraction and thus 
the oscillation frequencies of the anti-phase modes are strongly reduced and converge 
towards the values of the corresponding in-phase modes. This is found for the V2~ 
and V2+ modes around d — 2 and for the SR~ and SR+ modes for d > 2.182. As 
a consequence of the layer decoupling, the LR~ anti-phase layer rotation becomes the 
energetically lowest of the anti-phase modes. This indicates that the primary mechanism 
of decoupling of the electron and hole layers is the interlayer rotation LR~. 

4. Ground states and single-particle spectrum of quantum bilayers 

In this section we present an extension of the classical results of section 2 to quantum 
bilayers. Here, in contrast to the classical simulations, the ground state kinetic energy 
does not vanish even in the limit of temperatures T — > resulting in a finite spatial 
extension of the particle orbitals on the scale of the whole A^-particlc cluster. Hence, 
fermionic quantum features such as exchange effects (Pauli exclusion principle) must be 
included. 

In order to treat the e-h bilayer system of equations (1-3) quantum mechanically, we 
introduce the dimensionless coupling parameter A of a harmonically confined quantum 
system which relates the characteristic Coulomb energy Ec = e'^/{A7rexo) to the 
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characteristic confinement energy Eq — tvjJo 

~E^^ 1 1 — = — ' (6) 



where Xq — yjfij (mujo) denotes the oscillator length and as — 'ineh'^ / (me'^) is the 
effective electron (hole) Bohr radius. Thus Hamiltonian (4) including the kinetic energy 
can be rewritten in dimensionless form 

^ i=i i<j yy^i ~ ^j) 

using the transformation {r — > r/xo, E — > E/E^, d* d/xo}. Note that r and d* are 
measured in units of Xq and thus explicitly depend on the confinement frequency t^o- 
The characteristic energies and length scales of the classical (section 2) and quantum 
system are related by 

§ = (AV2)^/^ ^ = (2A)V3, (9) 

£/Q Xo 

so that the layer separations used in the Hamiltonians (3) and (8), respectively, are 
relatedby d* = (2A)V3rf. 

In the hmit A — > 0, both, electrons and holes behave as an ideal trapped Fermi 
gas independent of the layer separation d*. For A — >^ oo, it is Xo/a^ ^ 1, and quantum 
effects vanish. Thus one recovers classical behavior and shell configuration changes which 
coincide with those in table 1. At finite A, however, intra- and interlayer interactions, 
together with the parabolic confinement, give rise to a complex quantum many-body 
problem, which is subject of the following investigation. In the considered quantum case, 
ground state properties depend on the two parameters d* and A. Therefore, the question 
arises if the additional degree of freedom will induce additional structural changes. To 
answer this question, we performed self-consistent Hartree-Fock (SCHF) calculations of 
two coupled electron and hole layers of zero thickness, which are discussed in the next 
two subsections. 



4.1. Second quantization formulation 

In order to derive mean-field type equations for the e-h bilayer, we rewrite the exact 
Hamiltonian (7,8) in the second-quantized form H — He-\- Hh — He-h-, where 

^eW= y"dV^t^,)(r)/io(r)4w(r) (10) 
+ \j /d'^d'^V'I(^)(r)Vi^(^)(f) ^=L=4(;,)(f) , 

He-h dVd¥^t(r) (r) _ ^ ^ 4(f) 4(r) , (11) 
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with hQ{r) = |(— + r^) denoting the single-particle energy. Further, ■0g|)j)(r) is the 
annihilation (creation) operator of spin-polarizcd electrons and holes at space point r 
which satisfy the fermionic anti-commutation relations [■(/'e(/j)(r), ■(/'^^.^^(f)]^ = 6{r — r) 

and [■?/'g|]^j(r), ■?/'g|]^^(f)]^ = where [A, = Ab + bA. In a Hartree-Fock approach 
[41], the four field operator products entering equations (10) and (11) are approximated 
by sums over double products '4'l(h)''Pe(h) weighted by the generalized electron (hole) 
density matrix /7e(h)(i", f) = {'4'l(h)i^)''i'e{h){'^))e{h), where the expectation value (ensemble 
average) is defined as {A)e{h) = Trpe(/i)^. More precisely, with r?, ^ G {e,h}, the 4- 
operator products are approximated as 

i^liv)^{v)M^)i;,{v) (12) 

~ + Pr,(r,r)V'J(f)V'^(f) + p^(f,f)ViJ(r) V'^(r) 

- ^v^ [pr,(r,f)ViJ(f)V'^(r) + p^(f,r)V'J(r)V'^(f) 

Here, the first two terms constitute the Hartree term whereas the last two denote the 
Fock (exchange) contribution. The Kronecker delta assures that there is no exchange 
between electrons and holes which is due to the different physical nature of electrons and 
holes (different energy bands). Inserting the approximate expression (12) into (10,11) 
allows for an effective one-particle description according to 

Heih) =|y"dWf^^(,)(r){/.o(r)5(r-r) + S«^)(r,r)}4w(i-) , (13) 

^e-. =// dVd^f Vi|(r) {Ef^Jr, f) + E«^,(r, f)} 4(r) , (14) 
with the Hartree-Fock (HF) self-energies 

(r, r) = A / dV ^M!^ ,(r _ r) - A ^^^ptll , (15) 

^"^-HH e)(r, f ) = A /" dV /"(^^^""^''^^ 5{v - f) . (16) 
e h{h e)y , ; J ^(r' - r)2 + d*^ ^ ^ ^ ^ 

For computational reasons it is convenient to introduce a basis representation for the 
electron (hole) field operators, 

V'e(i)W- E^l*Hr)6i;i),, ie {0,1,2,...}, (17) 

i 

where the one-particle orbitals or wave functions (pi{r) form an orthonormal complete set 
and ag|)j-| ■ is the annihilation (creation) operator of a particle on the level i. Applying the 
basis expansion (17) to the equations (13) and (14) leads to the matrix representation 
of the bilayer Hamiltonian (1) which will be given in the following section, cf. equations 
(18-20). 
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4.2. Self- consistent Hartree-Fock simulation technique 

In matrix representation, the mean-field Hamiltonian for the bilayer system 
corresponding to the initial equations (1-3) reads 

/.y^) =hl + /.-^C^-'*) - , (18) 

^-eiH-H) ^ ^ J2 (wl'^t'^ - Wt,:j'-'^) ) , (19) 
kl 

^e-M.-e)^^^^e-.(.-e)^Me)^ (20) 



kl 



with the single-particle (orbital) quantum numbers i and j {k and /), hf,^^ being the 



electron (hole) total energy, hP^ - the single-particle (kinetic and confinement) energy and 
e{h h) ^^e^ h{h e)^ intra (inter) layer interactions in mean-field approximation. 
Further, p^j'*^ = i ^e{h),j) denotes the zero-temperature density matrix of electrons 
and holes with respect to the one-particle basis (pi{r). In equation (19) both the Hartree 
and the Fock contribution appear, whereas in equation (20) only the Hartree term enters. 
The explicit expression for the single-electron (-hole) integral is 

hl = l J drV*(r)(-V^ + r>,(r), (21) 

and the two-electron (two-hole) and electron-hole integrals are given by 

e-eih-h) ^ f [^2^^2- ^li^) ^li^) ^ji^) ^22) 

J J v^(r-r)2 + «*2 ' 

^T--" = 1 1 d'rd^r ^^^^ /^■(^)V^^-WV^' (^) , (23) 
"''''^ JJ V(r - r)2 + rf*2 

where a* ^ is utilized to avoid the Coulomb singularity for r — f . A small parameter 
of a* < 0.01 has been found to show convergence for all quantities of interest. Details 
will be given elsewhere [42]. 

For numerical implementation of the SCHF procedure yielding the eigenfunc- 
tions ^^'•'^''(r) (Hartree-Fock orbitals) and eigenenergies ef '^'' (Hartree-Fock energies) of 
Hamiltonian (18), we have chosen the orthonormal Cartesian (2D) harmonic oscillator 
states 

e-(^'+!/')/2 

^"^A^) = ,„ , , = '^rnix) Hniy) , (24) 

V 2"*+" ml nl TT 

with single-particle quantum numbers i — {m,n), r = {x,y), the Hermite polynomials 
Tim{x) and (m + l)-fold degenerate energy eigenvalues e^^„ = m -\- n -\- 1 where 
m,n & {0, 1,2,.. .}. The Hartree-Fock orbitals, expanded in the form 

'^fnr) = E4^'^^^W' (25) 

j=0 

with coefficients qj^"* G M and respective energies e1^^\ are obtained by iteratively solving 
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the self-consistent Roothaan-Hall equations [40] 

nt-l 

E - 4'^ 4^ - , (26) 

fc=0 

at fixed dimension n;, x n;, (i = 0, 1, . . . n;, — 1) according to standard techniques, for 
details see e.g. [41] and references therein. The resulting electron (hole) density p'jt\{Y) 
corresponding to given values of d* and A is defined as 

N-l N-l rii 

Pfii^) - E = E E ^^3%A ^'(r) , (27) 

fc=0 fc=0 /=0 

where each Hartrcc-Fock orbital k is occupied by a single particle only. 

For the e-h bilayers with N < 10, we used Ub = 50 of the energetically lowest 
oscillator functions 93m,n(i") to expand the Hartree-Fock orbitals, for = 12 we took 
rib — 55 which was sufficient to obtain convergent results. Due to the electron-hole 
attraction the cluster size is reduced compared to that of a single layer Coulomb cluster. 
This favours the use of a moderate number of basis functions to ensure convergence. ^ 

4-3. Transition from the ideal Fermi gas towards the classical limit 

Aim of this part is to investigate the transition from a strongly degenerate quantum 
system, i.e. A = 0, to the classical hmit A — > oo. To give a reasonable estimate for the 
range at which the classical ground state results become valid, we consider a system with 
= 6 electrons and holes at an intermediate layer separation of d* — 1.0. Of special 
interest will be the central spot of the (1,5) configuration which can most directly be 
assigned to a classical particle position. 

In contrast to the classical results the Hartree-Fock (HF) calculations fully take into 
account the wave nature of electrons and holes. The quantum many-body effects are 
evident already at A = 0. In the classical case, the total energy in the ground state is zero 
(all particles sit in the bottom of the trap). In the quantum case, this is prevented by the 
Pauh principle. Orbital-resolved HF calculations as function of coupling parameter A 
are displayed in figure 3. Here, the right panel shows the A'"-particle density p^}!\{v) and 

the six populated single-particle orbitals 0i^it\(i") for moderate (Ai = 5.0), intermediate 
(A2 = 15.0) and strong (A3 = 35.0) couphng. The SCHF results reveal that, in particular 
for small values of A, obviously several orbitals contribute collectively to the different 
high-density spots which unambiguously determine the cluster configuration. 

Concerning the lowest orbital i = 0, with an increase of A, the overlap with the 
higher orbitals vanishes and the wave function becomes localized when A exceeds a 
value of 35. In contrast, in the investigated range of A < 40 the other particles remain, 
independently of the observed density modulation, delocalized as can be seen on the 

f Note that the additional center particle in the case of A?" = 19 strongly increases the cluster size so 
that essentially more basis functions (nf, ^ 90) are required to ensure convergence. As the problem 
determining the two-particle integrals, equations (22) and (23), scales with C(n^) a computation is 
limited by memory requirements. 
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Figure 3. Ground state of the = 6 cluster as function of interaction strength A for 
fixed layer separation d* = 1.0. Right: Accumulated iV-particle density p^l'']_(r), on 

top, and corresponding single-particle HF orbitals (j^i^li* x(^) ^'-'^ three different coupling 
parameters A. The different signs of the wave function (blue and orange) are separated 
by white areas of zero amplitude, whereas areas of maximum amplitude are black. Note 
that the six high-density spots of the A^-particle density do not necessarily correspond 
to the single particles themselves as the configuration appears as a superposition of 
all orbitals. Left: (top) Electron-hole interaction energy equation (28), of the 

center electron and hole states for different approximations and (bottom) HF energy 
of occupied levels e^'"'*'' as function of A. 



orbital pictures. The transition towards the hmit of strong correlations can be estimated 
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from the e-h-interaction energy 

et(A) = -// dVdV-|0to(r)r l'^to(r)r , (28) 

of the electron and hole in the lowest orbital. The upper diagram in the left panel 
of figure 3 displays the A-dependence for four different approximations. For the ideal 
system, A = 0, electron and hole are not bound and 6^*7°'' vanishes. The black solid 
line shows the interaction energy (28) obtained from the SCHF simulations which for 
A <^ 1 agrees with perturbation theory (PT) where a linear A-dcpcndcncc follows from 
substituting the ideal wave function (pofi{r), see equation (24), for ^^i'g (r) in equation 
(28). 

For A » 1 a semi-classical result can be derived. Starting from the classical ground 
state configuration (1,5) the outer particles, together with the confinement, create an 
effective potential for both center particles which can be harmonically approximated. 
The direct quantum mechanical solution of the harmonic problem provides a finite 
Gaussian electron (hole) extension of width a = (jg = ah- Hence, the e-h-interaction 
energy (28) of the inner particles can be computed in a semi-classical way using 
0|io*(^) ~ (ct/tt)^^'' e~°'^^^"'"^^^/^. In the strongly correlated regime, starting at A > 30, 
the semi-classical and SCHF solution coincide very well. 

However, in an intermediate coupling range, A ~ 15, the e-h interaction energy is 
reduced compared to the semi-classical solution which reflects the fact that the orbital 
i — substantially deviates from a Gaussian, cf. the five side maxima of the orbital 
i = for A2 = 15 in figure 3. With increase of A this Gaussian becomes more and more 
peaked describing the transition to the classical limit \(f)'^^Q{r)\'^ 6{y).'^ Despite the 
good agreement with the semi-classical approximation, in the whole investigated range 
of A < 40 the system is found to be essentially non-classical. This becomes evident 
by comparing with the pure classical result e~^Zh = —X/d* which neglects any finite 
particle extension. Concerning all populated HF orbitals the transition towards the 
classical limit with increasing A is shown in the lower left diagram of figure 3 in terms of 
the orbital energies e^^^^ As mentioned in section 4.2 the harmonically confined ideal 
Fermi gas (A = 0) is (m -|- l)-fold degenerate with m,n E {0, 1, . . .}. Around A > 15 
the energy of the outer particles converges towards a five-fold degenerate energy which 
is separated from the (lower) energy of the center particle. 

4.4- Quantum ground state configurations and structural transitions for N = 10 

Beside the higher numerical effort of a single SCHF computation compared to its 
classical analogue, a complete study of the ground states requires, in addition to d 
and N, the exploration of A as a third degree of freedom. To overcome this problem 
and to reduce the task, we limit our investigation to the two-shell clusters iV = 10 and 
= 12 which were found to exhibit rich ground state properties in the classical limit. 

+ In the mean-field Hamiltonian (13,14) the classical limit is obtained by replacing pe(^)(r,r') 
'5(r-r')Ef=r S{r-r,). 



Coulomb-dipole transition in mesoscopic electron-hole bilayers 



17 




Figure 4. Left: (A,d*)-phase diagram for the iV = 10 bilayer in HF approximation. 
The configuration (TZ, 7) means delocahzation of charges on the inner ring TZ. The black 
sohd hne indicates the classical ground state transition (3,7) — * (2,8) which occurs at 
dcr = 1.0116ro from left to right. The red arrow points out an inverse transition 
compared to the classical (2, 8) (3, 7) crossing. The two dashed lines indicate 
the path when changing loq at fixed layer separation d for a germanium (e = 16eoj 
"^e(/i) ~ 0.25me) quantum-well structure, see equation (29); uji — 926 GHz, uj2 = 9.26 
THz, UJ3 = 98 GHz. 

Right: Electron (hole) density p'^l^\{v) at characteristic points marked (a) to (h) in the 
phase diagram. The side length of the contour plots is 9a;o- The open circles mark the 
corresponding classical ground state particle positions. The rightmost column displays 
the corresponding angle-averaged radial density profiles for d* — 0.5 (4.0), dashed (red) 
line. 



The analysis was done by systematically scanning the phase diagram for fixed 
values of d* ranging from 0.1 . . . 10.0. For each of these d* values we start from the 
ideal system at A = and increment the coupling parameter stepwise by 5X = 0.05. 
The convergence of each step is ensured by an adaptive, precision controlled iteration 
number with up to 2500 iterations of the Roothaan-Hall equations (26) per increment 6X. 
The described procedure allows for a systematic investigation of the phase diagram by 
a gradual transition from the ideal Fermi to the strongly coupled system. To verify the 
results obtained, the ground states with respect to individual points in the phase diagram 
were recomputed by starting from a random distribution as well as by decreasing the 
temperature of an initial (high temperature) thermal distribution [42]. All procedures 
are found to yield the same HF orbitals (energies) and thus the same A^-particle densities 
and shell structures. 

The results for the = 10 cluster are presented in figure 4. The ground state 
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phase diagram can be divided into four domains (left panel of figure 4): 

(i) At small A a weakly correlated degenerate Fermi liquid is observed within each layer 
(blue area in the left figure). The observed electron (hole) density is rotationally 
symmetric and exhibits non-monotonic radial modulations of an (nearly) ideal 
trapped Fermi gas. The proper density distributions for = 0.5 and — 4.0 
are shown in figures (a) and (e), respectively. 

(ii) At higher A two shells separate, see points (b), (f) within the red area in the phase 
diagram and the corresponding density profile (j). While on the inner ring the 
electron (hole) density is still isotropic, the density on the outer shell becomes 
angle-modulated and reveals seven high-density spots. The integrated position 
probability density on the inner and outer shells is close to 3 and 7, respectively. 
The configuration will be referred to in the following as (7^, 7) as on the inner ring 
TZ no localized density peaks in Prfl^]^(r) are present. Hence the nomenclature does 
not indicate the particle numbers, but the number of distinct density peaks, as the 
particle orbitals are delocalized over the entire cluster, see discussion in section 4.3 
and 4.6. 

(iii) Further increase of the coupling parameter leads to more pronounced (concentric) 
shells. In particular, the inner radial density decreases which is accompanied by 
the formation of angular density modulation, see figures (c) and (g). The shell 
configuration is found to be (3,7). 

(iv) At a certain value \cr{d*), the bilaycr system jumps from the (3, 7) into the (2, 8) 
shell configuration (green area in the phase diagram), see figures (g) — >• (h). 

The general behaviour of (i) - (iii) is independent of the layer separation d*. The 
localized (3,7) configuration (iii) emerges in two steps by rotational symmetry breaking 
from the Fermi liquid (i) maintaining a higher density on the inner than on the outer ring. 
However, an increase of d* beyond unity leads, by weakening of the interlaycr attraction, 
to a repulsive intralaycr and thus Coulomb-dominated coupling. Consequently, the 
cluster size increases, compare the density plots of figure 4 (a) vs. (e), (b) vs. (f), etc. . 
Moreover, for a fixed A » 1, the dipole-to-Coulomb transition towards the strongly 
correlated Coulomb regime induces the (2,8) shell configuration [28] which is observed 
when d* is increased from 0.5 to 4.0, see figure (d) vs. (h). This transition reduces the 
inner-shell density, see figure 4 right (red vs. dashed lines). 

Further, at a fixed d* > 2.0 an increase of A leads to a purely coupling induced 
configuration change (3,7) (2,8), see details in section 4.6. For d* = 10, both layers are 
already weakly coupled and become completely decoupled when d* is further increased. 
Consequently, the critical (blue, red and green) curves in the phase diagram converge 
towards horizontal lines. Note, that d* is measured in units of xq and thus depends on 
the confinement frequency cuq. This implies for an experimental setup, e.g. a double 
quantum-well heterostructure with fixed physical layer separation d, that one traces 
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hyperbolas of the form 

when changing the trap frequency ujq, see the dashed hnes in the phase diagram of 
figure 4. The larger the physical layer separation d (or effective particle mass fn*^^^-^), 
the more the hyperbola shifts to larger values of d* . Interestingly, e.g., for a germanium 
based quantum well, at fixed layer separation d = 1375A, the ground state structure of 
the quantum bilayer can be externally controlled by change of ujq only. 

A comparison of the classical particle positions (open circles in figures (a)-(h)), 
according to equations (9), with the shells and high-density spots of the HF calculations 
plotted in figure 4 reveals a good agreement. Larger cluster sizes compared to the 
classical case for small A are explained by repulsive fermionic exchange interactions. 
Further, the bold black line in the phase diagram indicates the classical transition from 
(3,7) to (2,8) which occurs at dcr — l.OllGro when crossing the line from left to right. 
It is found that the classical line gives a reasonable estimate also for the transition 
in the quantum bilayer system. Hence the trend, found in section 2 for the classical 
bilayer system, of center density reduction with increasing d also holds in the case of 
a strongly correlated quantum system, where the orbitals extend over several classical 
particle positions. In the classical limit, i.e. at very large A (outside of figure 4), the 
configuration boundary (3, 7) <-> (2, 8) (green curve) and the classical result (black 
curve) converge. Nevertheless, for intermediate values of A the red arrow indicates a 
remarkable point in the phase diagram where the structural transition in the classical 
and quantum bilayer proceeds in opposite direction. The single-particle orbitals for this 
transition will be analyzed in section 4.6. 

Further, an unusual (2,8) configuration is shown for A4 = 12.0 and rfg = 4.0 in 
figure 4(h), where the particle arrangement differs from the classical system. Such a 
configuration was also found in reference [28] for a classical single layer system with 
1/r" pair interaction and a < 0.94. Thus, the anomalous configuration underlines the 
effect of the Fermi repulsion in addition to the intralayer Coulomb interaction. However, 
an increase of A leads to a reduction of the Fermi effect and wave function overlap and 
a (2,8) configuration corresponding to the classical one is found. 

4-5. Quantum ground state configurations and structural transitions for N — 12 

In figure 5 we present the (A, d*)-phase diagram for N — 12 electrons and holes analogous 
to figure 4 for N — 10. At fixed (physical) layer separation d = lOOA one passes through 
four different domains of the phase diagram when the trap frequency is decreased from 
002 towards uji (see left panel of figure 5): 

(i) Analogously to the = 10 cluster at small A, a weakly correlated circular 
symmetric Fermi liquid exists within each layer, see point (a) in the blue area 
of the phase diagram. 
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Figure 5. (A,d*)-phase diagram showing the quantum shell structures found for 
the TV = 12 bilayer in HF approximation. The shown electron (hole) densities 
p^l'^^(r) corresponding to tuples (d) (d^Ai) = (0.5,5.9), (c) (d^,A2) = (0.75,3.95), 
(b) (dgjAs) = (1.25,2.35) and (a) (c?4,A4) = (2.0,1.45). The frequencies Wi^2 are as 
indicated in figure 4. The two black solid lines indicate the classical configuration 



transitions (3,9) (4,8) and (4, i 



(3,9) at Scr = 0.9528ro and d>cV = 0.3253ro, 



(2) 



respectively, from left to right. In the investigated range A < 15 these transitions were 
not observed in the quantum bilayer. The right two columns show the (radial) density 
of the four points (a) to (d) marked in the phase diagram for d = lOOA. 



(ii) A decrease of the trap frequency to point (b) is accompanied by a structural change 
to a 6-fold rotational cluster symmetry with an outer shape exhibiting hexagonal 
symmetry. This phase only establishes in the regime of a short-range in-layer 
potential, i.e. d* < 2. In the Coulomb case of weakly coupled layers this liquid-like 
state is not found. 

(iii) If the confinement strength is further reduced, see point (c), the cluster passes over 
to a 9-fold rotational symmetry. While in the cluster core a ring TZ of delocalized 
density is observed, the outer nine high-density spots are situated on a perfectly 
circular ring, which reproduces the symmetry of the external confinement potential. 

(iv) In the limit of small d* and A — > oo, see figure (d), where the in-layer interaction 
becomes extremely short-range, a commensurate closed packed structure with 
3— fold rotational symmetry similar to that known from classical dipole systems 
[28] is found. 

Consequently, during the coupling-induced transition from (i) to (iv) the cluster size 
decreases slightly as the effective in-layer interaction becomes short-ranged. In analogy 
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Figure 6. Left: Hartree-Fock energy eigenvalues ef-^^ corresponding to the spatial 
orbitals 0^^jt''^(r). (a) Bilayer density p^l'']^(r) and orbitals for = 10 at d* = 3.0 and 
Ai = 12.0, (b) A2 = 13.0. The black dots denote the orbitals which contribute most to 
the inner-shell high-density spots. While the inner shell of the (3,7) configuration is 
essentially build up from the 3 highest orbitals 7, 8 and 9, the inner shell of the (2,8) 
configuration is mostly formed from the orbitals 3 and 8 
Right: Single-particle orbitals 
of the wave function (blue and orange) are separated by white areas of zero amplitude, 
whereas areas of maximum amplitude are black. 



^jt\(r) for the cases (a) and (b). The different signs 



to = 10, the liquid-like state (i) as well as the (7^, 9) configuration (iii) are found for 
all values of d* . The additional configuration (ii), missing in the case of A^IO, is limited 
to a range of strong interlayer attraction. 

In contrast to = 10, in total two transitions as function of d were found in the 
classical = 12 system, cf. table 1. However, in the investigated quantum regime, 
A < 15, we observe no configuration changes corresponding to the classical transitions 
(3, 9) (4, 8), see black lines d^cr = 0.9528ro and d'^cr = 0.3253ro in figure 5 left. Hence 
the two ground state transitions (3, 9) (4, 8) and (4, 8) (3, 9) of type (A) and (B), 
introduced in section 2, are expected to occur outside of figure 5 in the (semi-) classical 
region only. 

4-6. Single-particle orbitals and single-particle spectrum 

In both previous subsections we discussed the phase diagram based on the A^-particle 
densities. In this part we pursue the question of how the single-particle spectrum evolves 
during the transition from (3,7) to (2,8) for the A^ = 10 cluster, see red arrow in figure 4 
left. At fixed d* = 3.0, the configurational transition occurs when changing the coupling 
parameter from Ai = 12 to A2 = 13. For this transition, the spatially resolved orbitals 
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A'"-particle density p^J*^l{r) are collected in figure 6 together with the 
corresponding one-particle HF spectra e^^^^ for both coupling parameters Ai and A2. 

As mentioned in section 4.4 the configuration change (2,8)ai ^ (3- T)\2 is reversed 
along the red arrow in figure 4 compared with the respective classical transition. Similar 
to the = 6 cluster discussed in figure 3 the HF orbitals generally extend over several 
classical particle positions. 

In situation (a), i.e. Ai = 12, the energetically highest orbitals i — 7,S and 9 contribute 
most to the inner-shell density showing three high-density spots. On the other hand, in 
(6), i.e. A2 = 13, the orbitals are completely rearranged with the two inner-shell density 
spots being now formed mainly from the orbitals 3 and 8, leading to embedded orbital 
energies e'^^'^^ and e^*^'*'* within the spectrum, cf. the black circles in the energy term 
schemes. In addition, all orbital energies of the (2, 8) configuration are enclosed in a 
narrower energy interval compared to (3, 7) whereas the energy spectra do not reveal 
any degeneracy. However, for (2,8) the spectrum separates into two parts of similar 
energetic substructure with orbitals energies e^^^^ to €4*-^^ and e^^'^^ to eg*-^^ respectively. 
Accompanying this fact, one clearly recognizes a change and an increase of the orbital 
symmetry when crossing over from the (3, 7) to the (2, 8) configuration. In contrast to (a) 
the rotational and specular (mirror) symmetry of 0i^jt\(r) with respect to perpendicular 
space axes in (6) is increased. Moreover, the structure of the nodes (white lines with 
zero amplitude in figure 6 right) of the HF orbitals changes, making the symmetry axes 
obvious. Particularly, inner and outer shell are clearly more separated by nodes in the 
(2, 8) configuration. 

5. Discussion and outlook 

In this paper we have considered ground state and dynamical properties of mesoscopic 
classical and quantum mass-symmetric electron-hole bilayers. In particular, we focused 

on the dependence of the properties on the layer separation d. The main effect is 
the gradual transition from systems with Coulomb interaction in the layers (at large 
d) to a system with short-range dipole interaction (at small d). Based on extensive 
classical molecular dynamics calculations we have shown that, with variation of d, several 
clusters show a sudden change of the ground state shell configuration, including several 
cases of re-entrant configuration changes which are related to symmetry properties. 
Furthermore, we have analyzed the classical normal modes of these bilayers and studied 
the d-dependence of the spectrum for A" = 19 as a representative example . 

A striking result is the energy jump of the inter-shell rotation mode frequency 0Ugj^+ 
by more than four orders of magnitude when the "magic" ground state configuration 
(1,6,12) is replaced by (1,7,11). This leads us to suggest a new possibility for external 
control of inter-shell rotation by exerting strain on the bilayer system (or alternatively 
by changing the trap frequency cuq by an external electric field [33]), i.e., a scheme which 
does not require changing particle number [43, 44]. Preparing a sample with d slightly 
above dcr, rapid compression initiates a ground state transition and thus allows to "turn 
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on" the intcr-shell rotation of composite dipoles — excitons. Combined with optical 
excitation this may have interesting apphcations manipulating coherent emission. 

In the second part of this paper we performed a quantum many-body calculation of 
the same system within the frame of a self-consistent Hartree-Fock approach. In the low- 
density limit where the particles are well localized the classical properties arc recovered. 
On the other hand, upon density increase and growing particle overlap quantum 
diffraction and exchange effects become important. This has significant consequences 
for the ground state phase diagram which is much richer than the classical one. There 
appear new structural phases which are characterized by charge localization on the outer 
shell coexisting with delocalization on the inner shell. Also, there exist parameter ranges 
where the classical and quantum systems show opposite shell configuration changes. The 
main advantage of the quantum many-body calculations is that they yield the complete 
single-particle energy spectrum and orbital-resolved ground states. We have shown that, 
even in the Wigner crystal phase where the density shows strong peaks, single peaks 
do not one-to-one correspond to single particles. On the contrary, in general, several 
orbitals contribute to a single density peak. 

We note that the present quantum results correspond only to the lowest level of 
many-body theory — the Hartree-Fock approximation. Thereby all pair interactions 
have been self-consistently included and direct and exchange terms are treated on the 
same footing. We have performed several comparisons with first-principle path integral 
Monte Carlo simulations which showed that the correct shell configurations arc observed. 
This lets us expect that the quantum results reported in this paper will not change 
qualitatively when better approximations are being considered. Naturally, the first 
improvement to be made is the inclusion of scattering effects on the level of the second 
Born approximation of nonequilibrium Green's functions theory, as it was done e.g. in 
references [45, 46, 47]. We presently develop these calculations which will be reported 
elsewhere. 
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